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Abstract. I discuss topics in the strong QCD dynamics that are of relevance to studies of 
the B mesons and to the spectroscopy of heavy hadrons. In particular the limitations from the 
hadronic dynamics for the determination of the weak-interaction parameters are discussed as 
well as some models for the newly observed apparently exotic resonances with hidden charm. 



1. Introduction 

The studies of the parameters of the Standard Model and the searches for a 'New Physics' in the 
properties of B mesons inevitably run into the usual problem that we formulate our theories in 
terms of quarks and gluons, while the experimental data are collected in terms of hadrons. Thus 
an interpretation of the wealth of data accumulated at the B factories within the underlying 
theory crucially depends on our ability to cope with the strong dynamics and to 'translate' 
between the quark-gluon and hadronic 'languages'. In this lecture I discuss some known results 
in such a 'translation' and the existing uncertainties, which still limit the precision of extracting 
the theoretical parameters. Some of these uncertainties can be somewhat relaxed with a help 
from new data, however some still appear to be beyond the current understanding of the strong 
dynamics. In what follows I discuss the current theoretical uncertainties in the precision of 
determination of the weak-interaction parameters from the data and, in some cases, the ways of 
reducing such uncertainties by further phenomenological studies of the yet unknown theoretical 
parameters resulting from the strong dynamics. 

The studies at the B factories thus far have not produced any credible indications of a 'New 
Physics'. However the unprecedented amount of high precision data accumulated in the process 
may be a treasure trove for learning new aspects of the QCD dynamics of quarks and leptons. In 
particular, one unexpected spin-off from the experiments at the B factories is the observation of 
a multitude of heavy resonances, both with open and hidden charm, displaying highly unusual 
properties, apparently indicating that some of those are the much expected exotic states, i.e. 
their structure goes beyond the standard quark-antiquark template for mesons. At present it 
looks like the new resonances fall into different categories according to their different internal 
structure. I discuss here two such types: the apparently 'molecular' state X(3872) and a class 
of resonances (Y, Z) which decay into specific charmonium states and pions. I describe in some 
detail the model, where these latter states are that of a 'hadro-charmonium', i.e. they can be 
considered as a compact charmonium state embedded in excited light-matter. 



2. Weak interaction parameters 

The weak decays of b hadrons provide us with the access to measurements of the entries in the 
Cabibbo - Kobayashi - Maskawa mixing matrix, describing the weak interaction mixing of the 
third quark generation. 

2.1. \V cb \ 

2.1.1. V c b from inclusive decay rate of B — > X c iv. The most straightforward technique for 
determining the mixing parameter \V c b\ is to apply the essentially parton picture to the rate 
of inclusive semileptonic decay of a b hadron, i.e. to identify the decay rate of a hadron, 
e.g. B — > X c £u, with the quark decay rate due to the process b — > civ, modulo small 
nonperturbative terms that take into account the 'Fermi motion' of the b quark and also the 
b quark spin interaction with the chromomagnetic field that the heavy quark 'sees' inside the 
hadron. Those nonperturbative terms are suppressed by the parameter ^qcd/ 171 ! an< ^ can 
be treated [JJ within an application of the Operator Product Expansion to the inclusive decay 
rates [2, 3j. Certainly, such approach assumes that the b quark mass is large enough for the 
uncontrollable low-energy non-perturbative effects to die out and thus for the quark-hadron 
duality to set in. It is not known at present to which accuracy this assumption works, however 
to some extent it is tested in the treatment of the lifetimes of b hadrons and their differences, 
as will be discussed further in this lecture. One trivial cross-check of such approach is to verify 
its prediction that in absolute terms the inclusive semileptonic decay rates B — > X c £v should be 
the same. The data agree with this expectation very well: according to the Tables [3], one can 
deduce T(B + -» £+vX) = 67.1 ± 1.8 ns" 1 T(B° -> £+uX) = 67.5 ± 1.8 ns" 1 . The charmless 
decay b — > ulv contributes only about 2% to each of these rates, so that this small contribution 
can be treated in either way without introducing an uncertainty exceeding the experimental 
error. The theoretical expression for the decay rate thus has the generic form (clearly related to 
the familiar muon decay expression) 

T(B — ► £vX c ) = — - — - — 5 — - / ( — - ) X [1 + (pert, correct) + (non — pert, correct)] , (1) 
1927r J Km),/ 

with the well known kinematical function 

f(x) = 1 - 8x 2 + 8x 6 -x s - 24x 4 lnx . (2) 

The explicit form of the corrections depends on the scheme used in a definition of the quark 
masses mj and m c [5j [6j [8] . The latter mass parameters and the uncertainty in them in fact 
dominate the problem of extracting \V c b\ from the data. One approach |9j to this problem is to 
fit these kinematical quantities from the data on the spectra in the semileptonic B decays, e.g. 
from the moments of the charged lepton energy Ef M n = f E™ ^f 1 dEg (in this classification 
r s z = Mo). The difficulty of this approach is in the nonuniform over the spectrum experimental 
sensitivity and in that higher moments are generally less certain theoretically due to larger 
corrections and also due to enhanced potential importance of the unknown deviations from the 
quark-hadron duality. A typical fit to the available data along these lines yields (see e.g. in |10| ) 
\Vd,\ ps (41.7 ±0.8) x 10" 3 , where the error includes the experimental errors and the error of the 
fit but does not include any theoretical uncertainty. 

2.1.2. \Vcb\ from exclusive decays B — > D^lv at zero recoil. An alternative way of extracting 
the mixing parameter | V c b \ is provided by the prediction from the heavy quark symmetry for the 
form factors of the decays B — > D*£u and B — > Dlv at zero recoil [TT], i.e. when the final charmed 
meson is produced at rest. Namely at the kinematical point where w = (pb ■ Pd) / '{MbMd) = 1 
the axial form factor jP: 

(D*(e)\(c l5 ^b)\B) = e,F(l) (3) 



and the vector form factor G: 



(D\(cr,ob)\B) =G(1) 



(4) 



are both close to one with only perturbative correction and corrections to the heavy quark limit: 
F(1),G(1) = 1+ pert.corr. +0(m" 2 , m^ 2 ). The averages |10| of the experimental results for 
the measurable products of these form factors with \V c b\ are F(l) \V c b\ = (35.97 ± 0.53) x 10~ 3 
and G(l) \V c b\ = (42.3 ± 4.5) x 10 -3 . Thus the real issue is the theoretical estimate of how far 
the form factors deviate from one as a result of the corrections. The axial form factor F(l) 
is estimated to be about 0.92 both analytically [5 J and numerically (on a lattice) [12] with an 
error believed to be of few percent, where 'few' typically implies 3 to 6. This estimate results in 
\V c b\ = [39.1 ± 0.6 exp ± (1.1 ~r~ 2.2)th] x 10~ 3 . The vector form factor is estimated numerically [13] 
as G(l) « 1.07, which results in |V cfe | = [39.1 ± 4.2 0xp ± 0.9i at ] x 10 3 . 

Although the two exclusive determinations of \V c b\ are in a perfect agreement with one another 
and also agree, within the uncertainties, with the inclusive determination, the estimates for the 
form factors may raise some reservations. The reason is that additionally to the heavy quark 
symmetry the closeness of G(l) (but not of F(l)) to one is also protected by the Ademollo-Gatto 
theoremp3] in the limit m c — > nib- In the real world the masses of the charm and bottom quarks 
are of course different, however the relevant parameter breaking the symmetry limit in G(l) is 
£ 2 = [(m& — m c )/(mb + m c )] 2 « 0.25 [11], and one could expect that some 'remnants' of the 
symmetry suppress the difference \G(1) — lj in comparison with the 'unprotected' |1 — F(l)\. 
Clearly, the currently used estimates of the deviations of the form factors from one show no traces 
of such a suppression. It can be noticed that if the actual values of the form factors indeed satisfy 
the Ademollo-Gatto requirement literally: G(l) — 1 ~ [1 — .F(l)], all three determinations of 
\V c b\ would naturally 'click' into place, noticeably better than they are currently believed to. 
Namely, with t; 2 = 0.25 one finds 
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in a perfect agreement with the determination from the inclusive semileptonic decay. The data 
then also imply that F(l) « 0.88 and G(l) « 1.03. 

It can be also mentioned that an agreement between all three determinations of \V c b\ can be 
viewed as a test of the V — A structure of the b — > c weak current. In particular the limit on 
the ratio of the right- and left-chiral amplitudes \V + -A|/|V — A\ is about 0.15 [15]. 



2.2. The semileptonic branching fraction B s i{B). 

Although it is not directly related to the determination of the CKM parameters, the low 
measured overall semileptonic branching fraction B s i{B) has been challenging to the theoretical 
understanding of the weak decays of the B mesons [16]. The nonperturbative effects (in the 
hadronic decay) are smaller for the B° meson. Thus the experimental number appropriate for 
a use as a reference is B(B° -> l + vX) = (10.33 ± 0.28)%. On the theoretical side B 8l (B°) is 
determined dominantly by the 'parton' decay branching ratio: 



. r(6 — ► ceu) 

sl ( -» ceu) = ^— - - r(fc - - ^ - r(ft - c _ d(s)) - — - c _ s(d)) - r(fc - - - ^ 

(6) 

with the main factor determining the ratio being the ratio of the quark and the semileptonic 
decay rates of the b quark: 

T(b -> cud(s)) _ a s a 2 / 2 15 \ 

3T(b^cev) " 1 + 2^ + ^ 4 [ L +Y L + 1 - 3 + h8 ) + --- ' (7) 



where L = \n(mw /m&) ~ 2.8. The calculation[17J of NNLO terms in the 0(a 2 ) radiative 
correction has increased the predicted relative hadronic decay rate and thereby lowered the 
prediction for the overall semileptonic branching fraction from 11.5% to 11%. As a result, the 
remaining gap between the theory and the data does not look as baffling [16] as it used to. 



2.3. Lifetime differences. 

The study of of the differences between inclusive decay rates of different b hadrons provides 
an insight into the applicability of using the heavy quark limit and the operator product 
expansion (OPE) in the Minkowski domain, i.e. the validity of the presumed quark-hadron 
duality. Within the OPE based approach the inclusive decay rate is expressed by means of the 
effective Lagrangian, which in fact is the absorptive part of the correlator of two (appropriate 
parts of the) weak interaction Lagrangians: 



"eff 



21m 



, j d A xe iqx T{L w {x),L w {<d)} 



(8) 



so that the decay rate of a heavy hadron Hq is calculated as 

V H = (H Q \L eff \H Q ) , (9) 

at q 2 = toq. For a heavy quark one can use the OPE in the inverse powers of the heavy quark 
mass niQ j2[ [3] : 

T _ r (°) i r ( 2 ) i r ( 3 ) - 

L eff ~ L eff + L eff + L eff ~ 

const • (QQ) + c 2 m% (q{3 ■ B)Q) + £ cf m 2 Q fer i% )(Q^Q) + . . . , (10) 

i 

where the first term describes the perturbative 'parton' decay, the second term introduces a 
correction due to the local gluonic field inside the hadron and the third term describes the 
effect of the spectator quark in the decay rate. The three first terms of the expansion in nig 1 
correspond to the unitary cuts of the diagrams shown in Figure [TJ 




Figure 1. The three types of graphs whose unitary cuts generate the first three terms of the 
heavy quark expansion for the inclusive weak decay rates. 



The difference in the inclusive decay rates of the b hadrons with different light-quark flavor 
arises from the third term of the expansion (110p . In particular one finds (3J [18] that the lifetimes 



of the B^ and B s should be the same within less than one percent, while the lifetime of the B A 
meson should be measurably longer: 



r(B±)-r(s°) = \v c 



d>\ 



|2 G p m b f B \,„n „r, „ _ 1 



8tt 



(Cl -C 2 _)B + -{C 2 + + C 2 _)B 



-0.030 



Sb \ 



2 



220MeVj PS 
(11) 

where fs is the B meson annihilation constant, the coefficients C+ and C_ are the standard 
non-leptonic weak interaction renormalization coefficients: 

C_ = C+ 2 = (a s (m b )/a s (m w )) 12 / 23 (12) 

and finally B and B are the 'bag' constants defined for a generic pseudoscalar heavy meson P 
as 

(PmKQ^il ~ 76)9) (97m(1 " 7 5 )Q)I*W = \f 2 P M P B 

(PQ<i\(Qj^ ~ 7s)Q) (97^(1 " 7?)g)|iW = \fpM P B . (13) 

In the so-called factorization approximation (where the four-quark operators are replaced by 
the products of two-quark with the vacuum insertion), both bag constants are equal to one: 
B = B = 1, in which approximation the numerical value of the decay rate difference in Eq. tjlip 
is estimated. 

Experimentally measured differences of the total decay rates are: T(B S ) — T(B^) = 0.027 ± 
0.013 ps~ l and T(B ± ) — T(B°) = —0.043 ±0.006 ps^ 1 , both in a reasonable agreement with the 
theoretical expectation, given the current errors. 

It can be also mentioned that the relative difference of the decay rates in the B hadrons is 
quite small, while similar differences for the charmed mesons are quite conspicuous. In particular 
by scaling the effect to from the B to D mesons, one estimates that the difference in the lifetimes 
between the charged and the neutral D mesons should be of order one: 

^m*4&^m„ om , 

T(D) m 2 c f% Y[B) 
which is very well known to be the case indeed. 



2.4-. \V u b\ and nonfactorizable terms in inclusive decay rates 

The determination of the small mixing parameter \V u b\ is complcated by the overwhelming 
background in the semileptonic decays arising from the dominant b — > c transition. In order 
to separate the charmless semileptonic decays B — > X u lu one has to impose kinematical cuts 
that eliminate the contribution of the heavier charmed final states X c . Such cuts however also 
leave only a minor fraction of the charmless decays as well. If, for example, the kinematical 
constrain is imposed on the value of q 2 of the lepton pairpl], q 2 > (M B -M D ) 2 , one is left with 
only about 20% of the total rate of the decays B — > X u lv. The theoretical calculation of such 
fractional decay rate into a restricted part of the phase space however suffers from uncertaintees 
which substantially limit the accuracy of the extraction of the mixing parameter \V u f,\ from the 
data[20j. In particular, within the OPE based approach the nonperturbative contribution of 
the third term in the effective Lagrangian (jlOh is formally concentrated in the 'useful' part of 
the phase space. (One can also formulate this behavior by recognizing that the effective mass 
parameter for the expansion of the decay rate restricted by q 2 > q 2 , is in fact /i = (m 2 — q®)/ (2mb) 
which becomes less than the mass of the charmed quark if q 2 , = (M B — M B ) 2 .) 



The nonperturbative correction to the semileptonic decay due to the b — > u current can be 
written as 



12 7T 

where Bi is still another bag constant defined as 

(B\b(l - 75 ) M )(u(l + 75 )6)|B) = (k^ B2 (is) 

and which also is equal to one in the heavy quark factorization limit. The actual values of the 
bag constants are currently not known. In the limit of large number of colors N the deviations 
from the factorization limit are expected to be suppressed by 1/N 2 , which puts the expected 
difference between the bag constants at O(0. 1) in the real world. Clearly, an 0(3%) uncertainty 
thus expected on the basis of Eq. (jl5|> is noticeable in comparison with the 20% usable rate in 
the determination of \V u b\- Thus an improvement in the precision of this parameter requires 
understanding the nonfactorizable terms. 

To some extent the nonfactorizable contributions to the decay rates can be probed by the 
experimental data. One quantity explained by these contributions is the lifetime difference 
between the D s and D° mesons. The observed approximately 25% diffrence agrees well with the 
understanding that the nonfactorizable terms should be of the order of 0.1 as compared to the 
factorizable ones (e.g. the difference in the lifetimes between the D and D°). A more direct 
access to the combination of the bag constants involved in Eq. (jl5p is provided by the difference 
in the total semileptonic decay rates of the D s and D° mesons. Namely one can find that 

r, i(D ») - r sim = o.08„- (^) 2 (^) 2 (g) , (17) 

where SB S is the SU(3) nonsinglet part of the difference {B2 — B) similar to the one in Eq. (|15p . 
Clearly, if 5B S ~ 0.1 the relative effect of this difference should be quite conspicuous given that 
the total semileptonic decay rate of D° is F s i(D°) ~ 0.165 ps^ 1 . For this reason a measurement 
of the difference in inclusive semileptonic decay rates between the D s and D° is of a great 
interest and may significantly reduce the theoretical uncertainty in determination of the mixing 
parameter \V u t>\, which currently is understood [10] to be in the range \V u b\ = (3 4- 5) x 10~ 3 . 



2.5. Some general considerations. 

The elements of the CKM matrix are currently known to a good extent. Of those that we not 
discussed here the parameter | Vt s | is essentially fully fixed by the unitariry of the mixing matrix: 
I Vi S I = \ Vcb\, and this value is in a perfect agreement with the measured frequency of the B s — B s 
oscillations, Am s = 17.8 ps -1 , while the parameter V t d is deduced from the Bj — Bd oscillation 
frequency, Am d = 0.51 ps -1 : \V t d\ ~ 8 x 10~ 3 , so that |Vtrf|/|Vtg| ~ 0.20. Combined with the 
measurement of sin 2/3 from the CP-odd time asymmetry in the B decays this in fact completes 
the test of the unitarity triangle and the measurement of its elements. 

It is often claimed that an improvement in the measurements of the CKM mixing parameters 
may probe a 'New Physics'. Although this is a logical possibility, it can be noted that the 
accuracy of the determination of the weak-interaction parameters as well other quantities related 
to quarks (e.g. masses) still remains limited by the theoretical understanding of the 'conversion' 
between the quark parameters and the properties of the observed hadrons. An illustrative 
example of this situation is provided by the original Cabibbo angle 9 C , i.e. the CKM parameter 
\V US \. The 'gold plated' mode for measuring it is the decay K — > ireu. The form factor 
/+(0) for this decay is fixed in the limit of the flavor SU(3) symmetry and is protected from the 



first order corrections due to SU(3) breaking by the Ademollo-Gatto theorem. The Kaon decay 
data are abundant and readily available, and the history of the measurements spans more than 
four decades. However I believe that the uncertainty in the understanding of the second order 
SU(3) breaking effects places a substantial intrinsic limit on the accuracy of the knowledge of 
sin# c . The 'natural' scale for the second order SU(3) violating effects is 0(10%). The Tables 
[4] list the value of \V US \ with a better than one percent error, based on the careful and detailed 
analysis of the value of /+(0) in Ref. [21]. The latter analysis however is based on the chiral 
perturbation theory and also uses a constituent quark model (on the light cone), which both 
have their limitations, especially if one aims at such a high accuracy. In particular some recent 
numerical (lattice) calculations do not confirm the results of the chiral analysis in |21j . (For a 
discussion see e.g. |4j.) Thus the current accuracy of the knowledge of the Cabibbo angle is not 
entirely clear. On the other hand it is also not clear what dramatically new conclusions could 
be drawn had we known sin# with a much better precision. At least the present author is not 
aware of any credible theoretical prediction of this parameter, which would be thus tested by 
the data. Neither it looks very likely that improving the precision in the Cabibbo angle and 
other mixing parameters would lead us directly to a 'New Physics'. After all — when was the 
last time a 'New Physics' has been found from the third (or even the second) significant digit 
in an otherwise known quantity, rather than from observing something qualitatively new? 

3. New hadronic dynamics 

Although extensive studies at the B factories have not brought any indications of a 'New Physics' 
they did bring something unexpected and extremely interesting. Namely, it appears that a 
whole realm of unusual states containing have quarks has been uncovered that may significantly 
extend the understanding of hadronic dynamics. Until quite recently the reliably known hadrons 
have all fit the standard quark model template: either quark-antiquark mesons or three quark 
baryons. Anything not fitting this scheme, generically referred to as 'exotics', was not seen 
with a convincing reliability, inspite of many searches. The situation has changed since the 
discovery[22j of a narrow (r < 2.3 MeV) resonance A(3872) produced in the decays B — > K X 
and decaying as A(3872) — > 7r + 7r~ J/ijj and with a mass in an extreme proximity to the D°D*° 
threshold: M X -M(D°D*°) = -0.6 ±0.6 MeV. By now it is pretty much clear that the A(3872) 
is a near-threshold singularity in the D°D*° S-wave channel with even C parity — either a 
shallow bound state [231 121 (23 [26] , or (more likely) a so-called virtual state [271 [281 [29] . Such 
'molecular' states made of pairs of heavy-light hadrons were expected since long ago \30\ 13 lj . 
Further studies at the B factories have uncovered new resonances with open and hidden charm, 
that exhibit exotic properties, although their internal dynamics is likely to be different from the 
'molecular' structure of A (3872). In the remaining part of this lecture I discuss some properties 
of the A (3872) and also of the so-called Y and Z resonances decaying into specific states of 
charmonium and pions, which may be states of an entirely new type, namely, of charmonium 
bound inside an excited light hadronic matter, 'hadro-charmonium' [32]. 

3.1. A(3872) 

The state A(3872) observed in B decays |22j [33] and in pp collisions at the Tevatron (34J, [35] . 
The co-existence of the decay A — > 7r + 7r~ J/ip (the discovery mode) on one hand and the decays 
X — > ir + ir~ir° J/ip and X — ► jJ/ip on the other implies that the isospin is badly broken in the 
A (3872). The C parity of the resonance has to be positive, and the bulk of data indicates that 
the quantum numbers of this state are J PC = 1 ++ . These properties can be explained if A is a 
C-even S-wave state of a charmed meson pair D°D*° + D°D*°. The threshold for a similar pair 
of charged mesons, D + D*~ + D~D* + is heavier by 5 ~ 8 MeV which is a large gap in the scale 
of the splitting of X from the D°D*° threshold. The presence in the wave function of the pair of 
charged mesons is thus heavily suppressed, which gives rise to A being a mixture of isospin and 
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Figure 2. The expected shape (in arbitrary units) of the virtual state peak in the yield of 7r + 7r~ J ftp 
(solid) and 7r + 7r~7r°J/ip (dashed) channels. The energy E (invariant mass) is measured from the 
D°D*° threshold. 



isospin 1 isotopic states. If X is dominantly a shallow bound state of the neutral mesons, those 
mesons should move inside the resonance at distances beyond the range of strong interaction, 
similarly to the proton and the neutron bound in a deutron. In this picture one should expect 
that the know decays of the D*°, or the D*° meson from the 'peripheral' component of the 
X(3872), D*° -»• jD°, D*° -> tt°D°, give rise to the decays of the X resonance X -» >yD D 
and X — ► 7r°D°D with a calculable rate and an interference pattern determined by the binding 
energy [24[ [36]. An experimental search for these decays produced an unexpected result: the 
peak in the jD D and 7t°D°D invariant mass spectra in the decays B — » XK turned out to 
be at or slightly above the threshold for D°D*° and was consistent with decays of unbound D*° 
mesons f37j [38]. This behavior suggests [28] that X(3872) is in fact a so-called virtual S-wave 
state in the C-even D°D*° + D°D*° channel, i.e. it corresponds to a large negative scattering 
length for the meson scattering near the threshold. In this case the scattering of the mesons 
near the threshold can be analyzed in both isotopic channels with 1 = and 1=1, and the 
interplay between the kinematics involving two thresholds for the D°D*° and D + D*~ states 
and the isospin symmetry of the underlying strong interaction produces a distinctive pattern 
of the dependence of the isospin structure of the scattering states on the energy E relative to 
the threshold [29|. This pattern can be tested by observing a different behavior of the threshold 
peaks in the X — > ir + ir~J/ip (I = 1) and the X — > ir + ir~ir° J/ip channels in the processes 
B — > XK as illustrated in Figure 2. 

3.2. Y,Z — hadro-charmonium? 

The charmonium resonances with mass below the open charm threshold (J /if},ip' ,Xcj, ■ ■ ■) ar e 
narrow due to suppression of the annihilation of a heavy quark pair into light hadrons. This 
suppression, known as the Okubo-Zweig-Iizuka (OZI) rule also finds its quantitative explanation 
in QCD [39]. Conversely, for higher charmonium resonances, i.e. those above the open charm 
threshold, the decay into pairs of charmed mesons should not be suppressed, and those resonances 
should be broad. The known charmonium states indeed generally follow this rule, e.g. ^(3770) 



is only 30 -f- 40 MeV above the DD threshold, and its width is about 25 MeV and is essentially 
all due to the decay into pairs of charmed mesons. However notable exceptions from such 
pattern have been recently found. The resonance Y(4260) with quantum numbers 1 has 
a large width of about 80 MeV, but non of its decays into charmed mesons have ever been 
seen. Rather the observed decay modes are Y — > mrJ/i/) and Y — » KKJ/ip [40], HU [42^ [43] . 
Subsequently more 1 peaks with similar behavior were observed in the range 4.32 - 4.36 MeV, 
Y(4.32 - 4.36) 011115] and near 4.66 MeV, F(4.66) @5], corresponding to the decay Y -> Tnrip' 
(but not J/tp). Finally, a manifest isotopic non-singlet Z(4430) has been reported[46j (but not 
yet independently confirmed [47] ) in B decays, corresponding to a resonance in the 7r ± ^ / channel, 
followed by an indication [38] of similar peaks Zi(4.05) and ^(4.25) decaying into tt XcX- 

The apparent existence of resonances that decay into a specific low-mass charmonium state 
and light hadrons invites an explanation [32J that these resonances in a sense already contain 
inside them that charmonium specific state. In other words, these high-mass resonances contain 
charmonium embeded in an excited light hadronic matter, so that in their observed decays the 
light degrees of freedom are 'shaken off'. Generally the idea of a binding of charmonium in 
a light matter is not novel — bound states of J /ip inside nuclei has been discussed for some 
time [491 150] . The novel feature of the structure of the possible new resonances, which we 
call hadro-charmonium, is that the ligh-matter 'host' for the charmonium is an excited (and 
unstable) state rather than a stable nucleus. The existence of such states can be understood in 
the limit of heavy quarkonium and of a highly excited light-matter resonance by considering a 
van der Waals type interaction between the quarkonium and the light matter [321 151] . 

The interaction of a compact charmonium state (call it generically ip) with long wave length 
gluonic field inside light hadronic matter can be described in terms of the multipole expansion 
in QCD[52] 153] 154] with the leading term being the El interaction with the chromo-electric field 
E a . The effective Hamiltonian arising in the second order in this interaction can then be written 
as 

H eS = -^E-Ef , (18) 

where Q W is the (diagonal) chromo-electric polarizability of the state ijj having the dimension 
of volume, and which can be expressed in terms of the Green's function Q of the heavy quark 
pair in a color octet state: 

« w = ^m a ngnC a m , (19) 
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with r being the relative position of the quark and the antiquark and £ a the difference of the 
color generators acting on them: £ a = t% — t% (a detailed discussion can be found e.g. in the 
review |32j). Finally, the QCD coupling g is included in the normalization of the field strength 
in Eq. (|18p . so that e.g. the gluon Lagrangian in this normalization takes the form — (l/4g 2 )F 2 . 

The diagonal chromo-electric polarizability is not yet known for either of the charmonium 
states, although it can be directly measured for the J /ip resonance [55j. What is known is the off- 
diagonal chromo-polarizability aft^ ) describing the strength of the amplitude of the transition 
ip(2S) —* ir + ir~ J/ip: \ct^'^\ ~ 2 GeV -3 [321 150], On general grounds, the diagonal chromo- 
polarizability for each of the J/tp, ip(2S) and Xc states should be real and positive, and for the 
former two the Schwartz inequality 

a W) a W>| a W)|2 (20) 

should hold. It can be expected however that each of the discussed diagonal parameters should 
exceed the off-diagonal one, and also that the chromo-polarizability for the ip(2S) and Xc states 
is larger than that of the J/ip due to their larger spatial size. 

The strength of the van der Waals type interaction between a charmonium state and a hadron 
(generically denoted here as X) can be evaluated using the effective Hamiltonian (I18j) and the 



well known expression for the conformal anomaly in QCD in the chiral limit: 

°" = "32^ F ^ Fa " V = 16^ {EtEt ~ W) ' (21) 

where B a is the chromo-magnetic field, and the (normalization) expression for the static average 
of the trace of the stress tensor 9^ over any state X in terms of its mass M x '■ 

(X\e£(q = 0)\X)=M x , (22) 

which is written here assuming nonrelativistic normalization for the state X: (X\X) = 1. 
Averaging the effective Hamiltonian (|18p over a hadron X made out of light quarks and/or 
gluons, one thus finds 

(X\H eS \X)<-^aM M x , (23) 



where the inequality arises from the fact that the average value of the full square operator BfBf 
over a physical hadron X has to be non-negative. The relation f|23|) shows an integral strength 
of the interaction. Namely, if the force between the charmonium {^) and the light hadronic 
matter inside the hadron X is described by a potential V(x), such that V goes to zero at large 
| a; |, this relation gives the bound for the integral 

V(x) d 3 x<-— a W M x . (24) 
9 

The value of the integral in Eq. (|24|) . although undoubtedly corresponding to an attraction, 
does not by itself automatically imply existence of a bound state, since it does not take into 
account the kinetic energy, which in a nonrelativistic treatment of the system is p 2 /2M ~ 
1/(R 2 M), with R being a characteristic size of the hadron X and M = M x / (M x + M$) the 
reduced mass in the system. The spatial integral of the kinetic energy is then parametrically of 
order R/M, and, given the relation in Eq. (l24l) . the condition for existence of a bound state can 
be written as 

a «) M^i>c, (25) 
R 

where C is a numerical constant (parametrically of order one) which depends on a model for 
the distribution of the interaction over the interior of the hadron X and thus on a more precise 
definition of the "characteristic size" R. Considering the condition (|25p for excited light-matter 
resonances, one can readily see that the appearance of a hadro-charmonium state, possibly at 
a higher excitation, depends on the behavior of the combination M x M/R with the excitation 
number of the resonance X. In particular, if the characteristic size R grows slower than the 
mass M x a binding of charmonium necessarily occurs in a sufficiently highly excited resonance. 
Such behavior takes place, e.g. in the once popular bag model[SB], where R oc M 1 / 3 . However 
in models which reproduce the approximately linear behavior of the Regge trajectories for the 
resonances, such as string model including its recently discussed [57J implementation in terms 
of AdS/QCD correspondence with 'linear' confinement, one effectively finds R oc M, and a 
more accurate estimate [58] proves that the binding of heavy quarkonium necessarily occurs 
at sufficiently high excitation of the light hadronic matter. Also one can argue [58] that in 
the limit of large mass itiq of the heavy quark the decay of such states into open heavy flavor 
channels, requiring a reconnection of the heavy- heavy and light-light binding, is suppressed as 

exp(- const. ^AqcWtoq). 

Obviously, the existing idealized theoretical picture can provide only qualitative guidance 
for the properties of the real-life hadro-quarkonium states. The most straightforward of these 
qualitative predictions are the following: 



• In addition to mesonic hadro-charmonium resonances there should exist bound states of 
baryo-charmonium, i.e. bound states of J /if), if)' , Xc with baryonic excitations at low baryon 
number. Such states should decay into e.g. a nucleon + J /if), possibly with additional pions. 

• There should exist resonances decaying into pion(s) plus other low-mass states of 
charmonium, r] c , %c- This prediction [32l[5T] is in a qualitative agreement with the apparent 
observation of Z\p [48] . 

• In addition to the observed decays into final states with the 'preferred' specific charmonium 
state, e.g. K(4260) — > ttttJ/i/j, the hadro-charmonium resonances should also decay, at 
a smaller rate, to non-preferred states, e.g. Y(4260) — > iririj)' due to a 'deformation' of 
charmonium in the bound state. 

• A similar to hadro-charmonium structure of resonances should exist for the bottomonium in 
the range around 11 - 11.5 GeV, most likely with the excited low-mass bottomonium states 
(due to a larger chromo-polarizability of the excited states). A possible hadro-bottomonium 
structure is in fact hinted by the data [59J on the inrT(nS) yield near 10.9 GeV. 
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